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Abstract 

The feature of the spectrum of primordial magnetic field is studied 
by using renormalization group analysis in magnetohydrodynamics. 
Taking account of the renormalized resistivity at the fixed point, we 
show that the scaling of the typical scale with time obeys L(t) ~ t 2 / 5 
for random initial condition. 



The magnetic fields observed in the various astrophysical scales is one 
of the important problems in cosmology @. The generation mechanism of 
magnetic field in proto-galaxy scale was already proposed by Kulsrud et 
al[|J. However, as this mechanism cannot explain the magnetic fields in 
intergalactic mediums and interclusters regions where matter is poor, it seems 
that the problem of the origin of magnetic fields still remains. This origin is 
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often attributed to the very early universe. Some attractive mechanisms in 
the course of cosmological phase transitions have been proposed In these 
mechanisms, the coherent length of the generated magnetic field is much 
smaller than the galaxy scale etc., because the magnetic field is generated 
by micro physical processes. Thus, one needs to investigate the process 
in which the small magnetic domains evolve to larger scales observed in the 
present universe. Since Kulsrud and Anderson || pointed out that the kinetic 
dynamo theory|| breaks down in interstellar mediums one cannot expect 
enough amplification of the magnetic field. Hence, we must re-analyze the 
evolution of MHD (Magnetohydrodynamics) from the starting point. 

Recently, Brandenburg et aljrj|| showed numerically that the inverse 
cascade occurs by using cascade model, and the occurrence of inverse cascade 
is also showed analytically by Olesen ||. In this paper, we reconsider Olesen's 
study from the viewpoint of the renormalization group. We show that the 
case considered in Ref. |J corresponds to our analysis at the fixed point. 
Then we discuss the evolution of the coherent length near the fixed point. 

Almost all of astrophysical systems have high magnetic Reynolds number 
because it is highly conductive in the typical astrophysical scales. Thus, 
MHD turbulence will appear in such system with fluid flow and we can use 
the statics of the turbulence to analize this system. We will apply the result 



obtained in Ref. ]I0 to the evolution of primordial magnetic field in our study. 

First of all, we review the renormalization group analysis of MHD. The 
basic equations of the incompressible MHD are given by 



and 



d t P + Q- VP = -Vp* + 7+AP + 7- AQ + f , (1) 



d t Q + P ■ VQ = -Vp* + 7+AQ + 7-AP + g, (2) 



where P = v + B, Q = v — B, p* = p + (1/2)B 2 , and j± = v ± rj. The f 
and g are the stirring forces and they are assumed to satisfy the statistical 
correlation, 

(fi{k)f 3 {q)) = 2k-yA (2n) 4 J ij (k)5 4 (k + q) (3) 
(fi(k) gj (q)) = 2k-yB«{2ityj l3 {K)5\k + q), (4) 
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and 



(9i{k) gj (q)) = 2k-yA (27i) 4 J tJ (k)5 4 (k + q), (5) 

where Jy(k) = Sij — kikj/\k 2 \ and A ,B Q are constants. The power of k 
should be decided by the initial condition. These forces are introduced so 
that the above system is equivalent to the original system in the inertial 
range( correspondence principle^llj). The solution is written as 

= /.w - i^ ( k)/ d5 ) , (6) 

where P,(k) = (gfl), />(*) = and J iJt (k) = %J it (k). G„(k) is 

the bare Green function written as 
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«o-'w=(- i ;: |k 7J k| _j;7 +|k| ,j. < 7 > 

Now, we divide the k— space into two parts, that is, large scale modes 
< k < Ae~ r and small scale modes Ae~ r < k < A, where A is the ultraviolet 
cut off. Then, we construct the effective equation for the large scale mode 
by integrating out the small scale modes. The result becomes 

G<-\k; r)P<(k) = f<(k) - zAo^(k) Jdq^ll &™ ) , (8) 

where subscript < means k e [0, Ae _r ], and G(k; r) is the dressed Green func- 
tion. Ao is the expansion parameter introduced technically. Now, y± depends 
on the renormalization parameter, r, and obey the following equation, 

where A d = l/(1207r 3 / 2 ), e = 1 + y, and F± are the functions of 7±, whose 
explicit expressions are not given here. From the above equation, one finds 
that the solution at the fixed point is y± oc e £r / 3 fT0[ . 



Compared the eq. (6) with eq. (7), one can see easily that the system is 
invariant under the following scale transformation, 

k -> k = ke r =:k£ uj -> Cj = ut a (10) 
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B(u£ a ,k£) -> B(k,u)) = B(k,u)r p v(w^°, k£) -»• v(M) = v(k,o;)r^(ll) 
v(t) ^v = u(i)i a - 2 ij(£) r? = r](£)r- 2 (12) 

f(k,^)^/(k,^) = f(k, w )r- /3 g (k, w ) - p(k,j>) = g(k,^r^. (13) 

From the requirement for the invariance of the stirring forces, one can obtain 
the relation, 

2(3 = 3a + (y + 3). (14) 
Now we can consider the evolution of the power spectrum E(k,t), 

(B 2 ) = J dkE{k,t) (15) 

From the scaling of (10)~(13), the equation 

£ 5+2a - 2/3 E(k, t) = E{£k, r a t) (16) 

holds. Defining the function ip(k,t) = k~ 5 ~ 2ot ~ 2 ^E{k,t), the above equation 
becomes 

ip(£k,r a t) = i>(k,t). (17) 
The solution is written as 

tp(k,t) = ^(k a t). (18) 
Thus, the spectrum is given by 

E{k,t) = k 2 ~ a - y ^{k a t) (19) 
As a = holds at the fixed point |TU[, we obtain 

E(k, t) = k {1 - 2y)/3 iP(k (5 - y)/3 t) (20) 

at the fixed point. The above expression is the same as one obtained by 
01esen||. Furthermore, as d t E(k,t) ~ E(k,t)/t ~ rj(k)k 2 E(k,t) approxi- 
mately at the fixed point, where r)(k) oc fc _<E//3 |I(|, the time evolution of the 
typical length obeys 

L(t) ~ kr x oc t 3 ^ 5 ~ y K (21) 
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Let us consider the case where the initial spectrum is white noise, y = 
—5/2, that is, E(k, 0) = k 2 ip(0). In this case, the typical correlation length 
becomes L(t) oc t 2//5 . If we assume that the magnetic field is generated in 
a cosmological phase transition at the temperature Tf (the cosmic time tf) 
and the coherent scale is the same as the bubble size, the present physical 
coherent scale is given by 

a{t f ) Hf' 

10 35 (A)(^ — r^GeV -1 ~ 1(T 2 x lOMpc, (22) 
U0- 4 ^100GeV^ * > \ J 

where lOMpc is the Silk damping scale of the magnetic fieldjl^] where the 
magnetic field with smaller scale is damped by the photon diffusion around 
the recombination epoch. is the ratio of the bubble size L(tf) to the 
horizon scale Hj 1 . For the electroweak phase transition, these typical values 
becomes /& ~ 10~ 4 and Tf ~ 100GeV[|l3j. This coherent length is 10 12 times 
as much as the comovingly developed scale; ~ 10AU. We note that this 
result differs from one obtained by numerical simulations in cascade model, 
L(t) oc t a25 0. However, since we used the renormalized resistivity which 
contains the non-linear effect fully, our estimation might be more useful. 

Next, to compare with the above result, we consider the case of Kraichnan 
spectrum E(k,0) oc k~ 3 ^ 2 , that is, L[t) oc t 4//3 , which is the one expected by 
the naive argument of Alfven effect in the inertial range ^M. In this case, the 



naive estimation shows that the coherent scale exceeds the present horizon 
scale as follows, 



L kra (t ) ~ l0 63 (-^-)(— L J—f GeV- 1 
v 07 Vio-4AiooGeW 

~ 10 28 x lOMpc > Hq 1 ~ 3000Mpc. (23) 



In this estimation, we need to study in the frame work of general relativity. 

From the above investigation, the correlation length can become enough 
size to explain the macroscopic scale. Although the Silk damping effect erases 
the structure smaller than the Silk scale, there is no drastic change of the 
quantitative result after taking account of the reasonable viscosity term||. 
Finally, we note that these results does not depend on the helicity of the 
magnetic field. 
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